Strain distributions in lattice-mismatched semiconductor core-shell nanowires 

Niels S0ndergaard 1 , Yuhui He 2 , Chun Fan 3 , Ruqi Han 2 , Thomas Guhr M , and H. Q. Xu 5 * 

1 Division of Mathematical Physics, LTH, Lund University, S-22100 Lund, Sweden 

2 Institute of Microelectronics, Peking University, Beijing 100871, China 

3 Computer Center of Peking University, Beijing 100871, China 

4 Fachbereich Physik, Universitat Duisburg- Essen, 47048 Duisburg, Germany and 

"'Division of Solid State Physics, Lund University, Box 118, S-22100 Lund, Sweden 



o 

c3 



O 






c3 



i 

C 

o 
o 



> 

in 
in 

in 

o 
o 

> 

• 1—1 

X 

5-1 



Abstract 

The authors study the elastic deformation field in 
lattice-mismatched core-shell nanowires with single and 
multiple shells. We consider infinite wires with a hexag- 
onal cross section under the assumption of translational 
symmetry. The strain distributions are found by min- 
imizing the elastic energy per unit cell using the finite 
element method. We find that the trace of the strain is 
discontinuous with a simple, almost piecewise variation 
between core and shell, whereas the individual compo- 
nents of the strain can exhibit complex variations. 



I. INTRODUCTION 

Nanowires have many applications. For example, 
nanowire-based biological sensors^, chemical detectors^, 
solar cells^™, LEDs^ field emitters^, transistors^, and 
electronic logic gates^- have been demonstrated. To fur- 
ther tailor the properties of nanowires, experimental- 
ists have grown nanowirc axial heterostructures, such as 
nanowire quantum dota^, and radial heterostructures, 
such as core-shelp2HHl anc j multi-shelf^ nanowires. In 
particular, a difference in lattice constant in heterostruc- 
tures leads to intrinsic strain, which for nanowires can 
be incorporated into the material much easier than for 
thin films due to the more effective strain relaxation at 
free surfaces^. This gives the opportunity for s train- 
engineering of the electronic and optical propertie d 17 * 18 !. 
The intrinsic elastic deformation field in heterostructured 
nanowires is, therefore, an important and necessary in- 
put for further investigations of the properties of these 
heterostructured nanowires. 

In this article, we discuss and calculate the strain 
field in lattice-mismatched core-shell and multi-shell 
nanowires. We find rich behaviour of the individual 
strain components, whereas their combination, the trace 
of the strain tensor (the volumetric strain) , shows much 
less variation. Our article is organized as follows. The 
model and theory is introduced in Sec. [H] The results 
and discussion are presented in Sec. |III| Finally Sec. [TV| 
contains the summary and conclusions. 
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II. THEORY AND METHOD 

In this work, we calculate the elastic deformation 
field in lattice-mismatched core-shell and double-shell 
nanowires (see Fig. fl] for schematics). As an approxima- 
tion, we take the nanowires as infinite with translational 
symmetry along the growth direction. We argue below 




FIG. 1: (a) Schematic of a nanowire with a single shell around 
a core, (b) Cross section of the core-shell nanowire with core 
radius n and shell radius r-z. (c) Cross section of a double- 
shell nanowire with core radius r\ and shell radii r-z and r$. 

(based on Saint- Venant's principle-^!) , that our descrip- 
tion of the infinite wires is also relevant for long finite 
wires. For simplicity, we neglect exterior forces acting on 
surface or bulk parts of the wires, i.e., we consider free 
nanowires. We restrict ourselves to nanowires with cu- 
bic lattice structure, but remark that generalizations to 
Wurzite or other structures are straightforward. Further- 
more, we consider only wires grown in the [lll]-direction 
in this work. For convenience, we define a primed coordi- 
nate system with respect to the growth direction, having 
basis vectors, 
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Consider, for example, a core-shell nanowire with un- 
deformed core and shell axial lengths of L c and L s and 
lattice constants of aS c ' and a^ s \ To allow for pseudo- 
morphic matching, we shall assume that both core and 
shell have the same number N of unit cells in the axial 



direction and, thus, L c must necessarily differ from L s 
We can write 



L^ 
a^ 






= C-N. 



(2) 



where C is some constant of proportionality. In the 
zincblende structure case for the [111] -direction, C would 
be v^- 

To match the lattice of the shell to the core we in- 
troduce the pseudomorphic initial strain field e*- ' in the 
shell. This choice of the pseudomorphic strain field ini- 
tially scales all shell lattice vectors to have the same 
length as in the core^l, 
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The total strain tensor is given by^l 

e= - (V(giu+ (V(g)u)*) 

For long b eams, an often used approximation is that 
of plane strai n 19 ' 22 '. In this approximation, the only non- 
zero strains are e' xx , e' and e' , whereas e' zx = e' = 
(planar sections remain flat) and e' zz = (no axial ex- 
tension). The theory has been generalized to stretched 
beams by assuming that the axial extension e' zz is non- 
vanishing but can be taken as a given fixed parameteif221. 
For the core-shell and the double-shell nanowires we con- 
sider in this work, e zz is not known a priori, but we 
expect it to be present due to the lattice mismatch. 
Therefore, we introduce an e' zz for each sub-domain of 
a nanowire and consider these strains as variables. In 
the simple case of a core-shell structure with two sub- 
domains (the core and the shell) , the matching effectively 

reduces the two axial strains e zz and e z z to a single vari- 
able a, 



p'W = J(i) - a 



a« 



— 1, with i = 



(6) 



In the material coordinate system, the potential elastic 
energy is 



U 



wdV= ~ J e-C-edV, (7) 



where w is the strain energy density, C — [Cij] is the 
matrix of elastic constants and e = [e^] are engineering 
strains 24 . In the nanowires with cubic lattice structure, 
the elastic constants are given by three independent con- 
stants taken by convention as C\\, C\ 2 and C44. 

From now on, we shall only work in the transformed 
coordinate system and drop all primes in the notation. In 
this coordinate system, we find, after some calculation, 
the energy density in Eq. Q as 

w = - (Di(e 2 -r . .,; L, y . ■ 



- ' el) + D 2 e\ + D 3 e 1 e 2 (8) 



-D 4 (ei +e 2 )e 3 + D 5 e\) , 



with the constants 

D 1 =: 

D 2 
D 3 
D 4 
D 5 



(Cn + C12 + 2C44) 



(9) 



= o(Ci 
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2C 12 + 4C 44 ) 
0C12 — 2C44) 
2C12 — 2C44) 



(C11 — C12 + 4C44) 



The potential energy allows us to find the deformation 
field by a variational principle. Minimizing the energy 
in Eq. Q leads to the equations for a collection of 
springs with generalized loads derived from the matching 
strairP2l. 

An important ingredient of this work is the introduc- 
tion of the axial extension as a variable. We shall there- 
fore briefly discuss the variational procedure for that de- 
gree of freedom. Now consider a core-shell wire. Axial 
variations obey 
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Varying the energy per unit cell with respect to the axial 
parameter a gives 
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dF, = F* 
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where we recognize that the axial stress is given in the 
form of 



a z 



dw 
de 3 
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Here we have used the assumption of translational invari- 
ance in the z-direction. Also, we note that in the second 
equality of Eq. (11), the integration in the axial direc- 
tion has been taken over the undeformed domains of the 
shell and core. In Eq. (11), a zz dS is the vertical force 



dF z on an area element dS. Hence, physically, Eq. (11) 



implies the vanishing of the total force in the axial di- 
rection. In the multi-shell nanowire case, the condition 



Eq. (11) generalizes to 



= V / a zz dS. 



(13) 



We remark that Eq. ( 13 1 should also apply to finite 
wires. Thus, using the principle of Saint- Venant-^, we 
expect our results to describe well the strain field in 
the middle sections of finite, but long, free wires. If a 
net total force and moment on the terminating ends are 
present, the corresponding condition Eq. (13) should be 
changed accordingly. 

In the numerics, we have used the line ar fi nite element 
method^ to minimize Eq. (I7| with Eq. (11) imposed for 
a free two-dimensional boundary. The planar degrees of 
freedom are seen to couple to the axial degree of freedom 
via, e.g., Eq. (ITT). 
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FIG. 3: Elastic strain distributions in a double-shell 
nanowire. The core with radius ri = 3.5 nm is made from 
GaP, the inner shell with radius ri — 9.2 nm is made from 
GaAs, and the outer shell with radius r$ = 13.9 nm is made 
from GaP. The plots show cross-section distributions of (a) 
e X x, (b) t xx , (c) e xy , and (d) tre. 
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FIG. 2: Elastic strain distributions in a single-shell nanowire. 
The core with radius n = 6.0 nm is made of GaAs and the 
shell with radius ri — 13.9 nm is made of GaP. The plots 
show cross-section distributions of (a) e xx , (b) e yy , (c) e xy , 
and (d) tr e. 



III. NUMERICAL RESULTS 

First, we consider a nanowire with a single GaP shell 
and a GaAs core with the same geometry as in the ex- 
periment reported in Ref. Q3J We choose radii of the 
inner core and shell to be r\ = 6.0 nm and T2 = 13.9 nm 
(see Fig. [TJfor the definitions of r\ and r^). We depict 
the strain components e xx (a), e yy (b), and e xy (c), and 
the trace tr e (d) in Fig. [2] (the strain component e zz is 
omitted). Figure [2] (a) shows that near a horizontal in- 
terface the exterior material is expanded. Figure [21 (b) 
shows that the exterior material is largely expanded in 
the regions near the left and right corners of the core 



of the hexagonal shape. These behaviours are consistent 
with the difference in the lattice constants. Figures [2] (a) 
and (b) also show that the two strain distribution pat- 
terns are symmetric with respect to the x- and the y-axis. 
The strain component e zz , which is not shown here, ex- 
hibits a simple step-like distribution profile with a con- 
stant axial contraction in the core and an axial elongation 
in the shell, again consistent with the lattice mismatch. 
The shear strain e xy in Fig. [2] (c) , which is important in, 
e.g., the piezo-electric interaction, displays a pattern as 
complicated as the other individual strain components, 
e xx and e yy . The distribution pattern of e xy is, how- 
ever, antisymmetric with respect to the x-axis and the 
y-axis. The measure of volume deformation, the dilata- 
tion SV/V — tre, is depicted in Fig. [2] (d). In general, 
the shell is expanded particularly near the boundary in 
the middle of the edges. This expansion extends typically 
all the way into the corners of the core. The interior of 
the core is seen in Fig. [2] (d) to be contracted almost 
constantly. 

We now turn to the case of multi-shell nanowires. For 
simplicity, we vary the setup only slightly by consider- 
ing a core with two shells, a double-shell nanowire. We 
choose the material of the inner shell to be made of GaAs 
and the core and outer shell are chosen to be made from 



GaP. Here, in consistence with the previous geometry, 
we choose the outer shell to be parallel with the inner 
core, see Fig. fl] (c). Still, we mention that fully parallel 
multi-shells are also of interest!^. The radii of the core, 
inner shell, and outer shell are taken to be r\ = 3.5 nm, 
r-2 = 9.2 nm and r% — 13.9 nm. Figure [3] shows the strain 
components e xx (a), e yy (b), and e xy (c), and the trace 
tr e (d) of the double-shell nanowire. (Here, again, the 
strain component e zz is omitted, since it shows a simple 
step-like distribution profile.) Since the lattice constant 
of the inner shell is larger than its surroundings, the in- 
ner shell is compressed whereas the core and outer shell 
regions are expanded (see plots for the strain components 
£xx, £yy, and the trace tre in Fig. pjj). In particular, the 
exterior shell near the interface undergoes considerable 
stretching as in the single-shell case. Also the same as in 
the single-shell nanowire case, the strain patterns of the 
components e xx and e yy are symmetric with respect to 
the x- and the y-axis, and the pattern of the trace tr e is 
hexagonal symmetric. Furthermore, the inner shell plays 
almost the same role for the tensile strains as the in- 
ner core in the single-shell case. The shear strain e xy in 
Fig. pi (c) shows similar symmetric characteristics as in 
the single-shell nanowire, i.e., the distribution pattern of 
the shear strain e xy is antisymmetric with respect to the 
x- and the y-axis, and exhibits similar peak- valley con- 
figurations in both inner shell and outer shell as in the 
shell region of the single-shell nanowire. 



IV. SUMMARY AND CONCLUSIONS 



We present a theoretical study of the strain field in 
lattice-mismatched core-shell nanowires with single and 
double shells. We derive a functional for the elastic en- 
ergy in a nanowire. The deformation field is found by 
minimization of the energy functional using the finite el- 
ement method. For the single-shell wire, the core is made 
of GaAs and the shell of GaP. For the double-shell wire, 
the inner shell is made of GaAs, whereas the core and 
outer shell are made of GaP. We find a large volumetric 
strain in various regions of a wire. A large compression 
appears in the core region of a single-shell nanowire or 
in the inner shell region of a double-shell nanowire. The 
large volumetric strain will influence carriers via the de- 
formation potential. Our numerics also shows great vari- 
ations in the individual components of the strain com- 
pared to the volumetric strain. This could be of particu- 
lar importance for the hole confinement in p-type wires. 
We note that the theory presented in this work can be 
extended to incorporate the finite strain components e xz 
and e yz which have been neglected in the plane strain ap- 
proximation. We leave the details of such a more general 
theory to forthcoming publications. 

The authors thank F. Boxberg for valuable discussions. 
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